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ABSTRACT 

This paper presents spectra in the 2 to 20 /im range of quiescent cloud material located in the 
IC 5146 cloud complex. The spectra were obtained with NASA's Infrared Telescope Facility (IRTF) 
SpeX instrument and the Spitzer Space Telescope's Infrared Spectrometer. We use these spectra to 
investigate dust and ice absorption features in pristine regions of the cloud that are unaltered by 
embedded stars. We find that the Il20-ice threshold extinction is 4.03 ± 0.05 mag. Once foreground 
extinction is taken into account, however, the threshold drops to 3.2 mag, equivalent to that found 
for the Taurus dark cloud, generally assumed to be the touchstone quiescent cloud against which all 
other dense cloud and embedded young stellar object observations are compared. Substructure in the 
trough of the silicate band for two sources is attributed to CII3OH and NII3 in the ices, present at 
the ^ 2% and ^ 5% levels, respectively, relative to Il20-ice. The correlation of the silicate feature 
with the E{J — K) color excess is found to follow a much shallower slope relative to lines of sight that 
probe diffuse clouds, supporting the previous results by Chiar et al. (2007). 

Subject headings: dust, extinction — infrared: ISM — ISM: lines and bands — ISM: molecules — 
ISM: individual (IC 5146) 



1. INTRODUCTION 

Ices are abundant in preplanetary material and their 
presence is likely to represent a vital resource for the ori- 
gin of planetary life (e.g., Mottl et al. 2007). Ices are 
widely observed in the dense molecular clouds that give 
birth to new stars (e.g., Gibb et al. 2004); they are ubiq- 
uitous in comets (e.g., Whipple 1950; Schulz et al. 2006), 
and have recently been detected on asteroidal bodies as 
well (Campins et al. 2010; Rivkin & Emery 2010). To de- 
termine the evolution of such materials, from icy coatings 
on submicron-sized interstellar grains to planetesimals in 
protoplanetary disks, is a major goal in the astronom- 
ical search for origins. A crucial first step is the for- 
mation of ice mixtures composed primarily of H2O, CO 
and CO2 in molecular clouds. The observed abundances 
of these molecules in the solid phase can be explained 
by simple grain surface reactions such as O + H — s> OH, 
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OH + H ^ H2O and CO + OH ^ CO2 + H, in combina- 
tion with direct freeze-out of CO from the gas phase. 

The physical conditions under which these reactions 
occur may vary from cloud to cloud, as a function of fac- 
tors such as the cloud mass, temperature, and the inten- 
sity of the local interstellar radiation field. Ices are quite 
widely observed in the spectra of young stellar objects 
(YSOs), but the study of ices in pristine prestellar dark 
clouds requires observations of background field stars in 
lines of sight that do not intercept circumstellar mate- 
rial around YSOs. Relatively few such observations are 
currently available. The dark cloud in Taurus is by far 
the best studied to date (e.g., Whittet et al. 1988, 2007; 
Chiar et al. 1995; Bergin et al 2005; Shenoy et al. 2008) 
and is often assumed to be prototypical, but this remains 
to be confirmed. Observations of ices toward a signif- 
icant sample of background stars are available for only 
two other molecular clouds: Serpens and p Oph (Eiroa & 
Hodapp 1989; Tanaka et al 1990; Chiar et al. 1994; Knez 
et al. 2005). Boogert et al. (2011) present a study of ices 
and extinction toward sixteen isolated cores as probed 
by background stars and find the threshold extinction 
to be in line with that found in Taurus. The observa- 
tions show similarities but also significant differences in 
ice properties and the conditions needed for ice to form. 
These differences may be parameterized in terms of the 
ice "extinction threshold," i.e., the minimum observed 
continuum extinction necessary for the detection of the 
ices (e.g., Whittet 2003), and in the relative abundances 
of H2O, CO and CO2 (Whittet et al. 2007, 2009). To es- 
tablish a reliable benchmark for studies of ice evolution 
in regions of active star formation it is vital to under- 
stand the initial conditions in quiescent clouds and the 
extent to which they vary from cloud to cloud. 

With this goal in mind, we present new observations 
of ices in the IC 5146 dark cloud complex. IC 5146 is 
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nearly 3° in extent, located in the direction of Cygnus 
(Elias 1978; Dobashi et al. 1992). The IC 5146 complex 
consists of several dark clouds as well as the well-known 
H II region, the Cocoon Nebula (Dobashi et al. 1992). 
The dark cloud called the Northern Streamer, a region 
of high extinction compared to the rest of the cloud, is 
located to the west of the Cocoon (Cambresy 1999; Lada 
et al. 1999; Harvey et al. 2008). It contains several pro- 
tostellar clumps, some of which are actively forming stars 
(Dobashi et al. 1992). These distinct regions are thought 
to be co-distant at 950 pc (Harvey et al. 2008). Star- 
formation in the H II region is loosely clustered and more 
efficient than that in the Northwest Streamer where most 
of the YSOs have formed in isolation. The field stars 
presented in this paper probe this relatively quiescent 
Northwest Streamer region. In fJDwe describe the com- 
plete data set for the IC 5146 cloud including the Spitzer 
(IRS, IRAC, MIPS), 2MASS (Skrutskic et al. 2006), and 
IRTF-SpeX observations. We used our IRTF-SpeX spec- 
tra to classify the spectral types of the background stars 
(i j2.1[) . thus allowing us to fit the continuum flux using 
a reddened photospheric model appropriate for the spec- 
tral type of the star fi i2.2|) . The resulting optical depth 
proflles and ice abundances are discussed in iJ3] The re- 
lationship between the near-infrared color excess and the 
silicate optical depth is discussed in fJH Finally, we sum- 
marize our results in Sj5j 

2. OBSERVATIONS AND DATA REDUCTION 

We used the Infrared Spectrometer (Houck et al. 2004) 
on the Spitzer Space Telescope (Werner et al. 2004) to 
obtain data in the 5 to 20 ^m region of 10 sources be- 
hind the IC 5146 dark cloud, with visual extinction (Ay) 
in the range from 4 to 27 mag. Source names, 2MASS 
IDs, spectral class (see §2.ip . extinction information, and 
AOR keys are listed in Table [1] The IRS pipeline version 
used was 15.3. The low resolution Short-Low and Long- 
Low modules (SL and LL; i? ~ 60 — 120) were reduced in 
a way that is customary for ground-based spectra, follow- 
ing the method described in Boogert et al. (2008). For 
the IC 5146 cloud, the background is relatively uniform 
across the array, so background subtraction was carried 
out by using the nodded pairs. A fixed width extraction 
was performed and the 1-dimensional spectra were then 
averaged. Subsequently the IC 5146 spectra were divided 
by the ratioed spectra of the standard stars HR 2194 (AO 
V; PID 1417; AOR keys 0013024512 and 0013024768 for 
SLl and SL2 respectively) and HR 6606 (G9 HI; PID 
1421; AOR key 0013731840 for LL) reduced in the same 
way in order to correct for wavelength-dependent slit 
losses. The standard stars were ratioed with the model 
spectra from Decin et al. (2007) appropriate for their 
spectral type. 

For seven of these sources, we also obtained 1.9 to 
4.0 /im data with the SpeX instrument (Rayner et al. 
2003) on NASA's Infrared Telescope Facility (IRTF). 
The spectra were obtained using a 0.8" slit with the 
cross-dispersed LXD mode of the SpeX instrument re- 
sulting in i? '--^ 937. The spectra were reduced using the 
SpeXtool software package (Vacca et al. 2003; Cushing 
et al. 2004). These data were used to classify the spectral 
types of the stars and measure the 3.0 /Ltm H2O ice ab- 
sorption feature. For the most heavily obscured source, 
Q21-1 we also obtained an M-band spectrum in order to 
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Figure 1. Spectral energy distributions for the IC 5146 field stars. 
2MASS and IRAC pfiotometry (diamonds), SpeX and IRS spectra 
(black lines) and the reddened photospheric model spectrum (green 
line) are shown. The spectra are plotted in order of decreasing 
visual extinction; Q21-1 {Av=23.9 mag), Q21-6 (Av=20.7 mag), 
Q22-1 (Av = 17.49 mag). 

measure the CO-ice profile at 4.67 /im. 

In order to flux calibrate the Spitzer-IRS and IRTF- 
SpeX spectra, we used available photometry from 
2MASS and IRAC (Harvey et al. 2008). The convolved 
flux at the effective wavelength of the photometric fil- 
ter was calculated using published responsivity curves 
for 2MASS (Cohen et al. 2003) and IRAC (Reach et al. 
2005). Following the calculation of the effective ffux over 
a particular photometric passband, a normalization fac- 
tor was calculated individually for the SpeX and IRS 
spectra. SpeX spectra were normalized to the 2MASS 
Ks band, and IRS spectra were normalized to either 
the IRAC3 or IRAC4 band. The flux calibrated spec- 
tra along with the photometry are shown in Fig. [TJ 

2.1. Spectral Classification 

In order to conflrm the background nature of our 
sources, we constructed color-color diagrams which sepa- 
rate out reddened field stars from embedded objects and 
stars with circumstellar material. Figure [5] shows the 
placement of the IC 5146 sources with respect to normal 
giant and dwarf stars (Bessell & Brett 1988)|3j along 
with a reddening vector. Reddening curves are deduced 
using the extinction law determined by Indebetouw et al. 
(2005) and are equivalent to Ay ^ 5 mag. It is appar- 

We chose not to transform these intrinsic colors to the 2MASS 
system because the transformation artificially displaces the "ori- 
gin" of the dwarf branch away from J — Kg = H — Kg = 0. 
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Table 1 

Field Stars located behind IC 5146 



2MASS ID 


Spitzer AOR label 


Spectral Type^ 




Ay from model match 


AOR Key 


21472204+4734410 


Q21-1 


K2 III Fe-0.5 


4.5+°-«^ 

^-0.07 


23.9 


14136064 




Q21-2 


na TTT Pp 1 


r, q + 0.05 

^•■^-0.00 


12.2 




21475842+4737164 


Q21-3 


[K2 III] 


2.1 


11.1 


10750720 


il45U774+47oll51 


Q21-4 


K6 111 OlN 0.5 


1 f.+0.08 
^•°-0.06 


8.5 


10750720 


21444787+4732574 


Q21-5 


K2 III Fe-0.5 


1 t:+0.07 
^•■^-0.07 


8.0 


14136064 


21461164+4734542 


Q21-6 


G8.5 Ilia Fe-0.5 


o n+0.05 

■'•^-o.oo 


20.7 


14136064 


21443293+4734569 


Q22-1 


KO III 


o q+0.04 

■^•■^-0.05 


17.49 


14136064 


21473989+4735485 


Q22-3 


[Kl III] 


0.7 


3.7 


10750976 


21473509+4737164 


Q23-1 


Kl Illb CN1.5 Cal 


1 4+0.07 
.^•^-0.04 


7.4 


14136064 


21472220+4738045 


Q23-2 


[K5 III] 


1.9 


10.1 


14136064 



^ Spectral typo based on mateh to standard star spectra in IRTF Spectral Library, with reddening applied. Square brackets indicate no IRTF data 
were available to classify the star and spectral type is based on fit of reddened Decin models to photometry and Spitzer IRS spectra. 



Intrinsic colors from Tokunaga 2000 
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Figure 1. Q21-2 (Av=12.2 mag), Q21-3 (Av=ll.l mag), Q23-2 
{Av=lO-l mag). 

ent from these diagrams that all the IC 5146 program 
stars are field giants reddened by dense cloud dust that 
displaces them from the red-giant branch. Young stellar 
objects and objects with circumstellar dust would fall 
near the bottom right in the J ~ Kg vs. [4.5] — Ks dia- 
gram. 

We used the 2.2 to 2.5 /im region to determine the 
spectral type of the stars for which we have IRTF-SpeX 
spectra (Fig. [3a|). In this spectral region, the strength 
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Figure 1. Q21-4 (Av=8.5 mag), Q21-5 (Av=8.0 mag), Q23-1 
{Av=7A mag). 

of the CO lines is a strong indicator of the spectral type 
(and luminosity class) of late-type stars. The ^^CO first 
overtone, spanning the 2.29-2.50 /im region, dominates 
the K-band spectra of the giant stars (Wallace & Hinkle 
1996, 1997; Heras et al. 2002) and its strength increases 
with decreasing temperature (i.e., later spectral types). 
The 2-0, 3-1, and 2-4 bands of the ^^G^^O isotope are 
also present in these giant stars due to the larger CO ab- 
sorption and decreased ^^C/^^C ratio (Hinkle & Wallace 
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Figure 1. Q22-3 {Av=^.7 mag) 

1996). Additionally, the Ca multiplet and the Fe line in 
the 2.261-2.267 /im region strengthen somewhat in spec- 
tral types later than early G. The Mg line near 2.28 /im 
does not show much variation with temperature (Wal- 
lace & Hinkle 2006), and is present at similar strengths 
in all our spectra. In order to classify our background 
stars, we compared our observed spectra to the stan- 
dard star spectra available in the IRTF Spectral Library 
which contains 57 GO through MIO giant-star spectra 
(Gushing, Rayner, & Vacca 2005; Rayner et al. 2009). 
The standard-star spectra were first reddened by apply- 
ing the Indebetouw et al. (2005) extinction law from 1.25 
to 8.0 /im, assuming constant extinction longward of 8 
/im. Wc assume Ak/Av — 0.09, appropriate for R=3.05 
(Whittet 2003). The visual extinction. Ay, was deduced 
based on E{J — K), assuming average intrinsic colors 
for late G through early M giants from Tokunaga et al. 
(2000) as a first best guess, and Ay = 5.3*E{J — K), ap- 
propriate for dense clouds (Whittet et al. 2001). Using 
the deduced extinction curve and taking advantage of the 
sensitivity of the GO lines to temperature, we determined 
the spectral type for each source (Fig. [3al Table 1). The 
uncertainty is no more than ±1 in subtype. Following 
the classification, E{J — K) was recalculated using the 
intrinsic colors for the appropriate spectral type; these 
values are given in Table 1. 

2.2. Determination of Continuum and Optical Depth 

Profiles 

The determination of an underlying stellar continuum 
is one of the biggest sources of uncertainty when deduc- 
ing ice and dust absorption profiles. Polynomial continua 
are somewhat arbitrary, but are often used when the true 
stellar continuum is not known and/or cannot be easily 
modeled. Since we have photometry and spectroscopy 
over a wide wavelength region and we are able to deter- 
mine the spectral types of our program stars with relative 
accuracy (^ N2.ip . we are able to fit reddened photospheric 
models from Decin & Eriksson (2007). These models 
cover the spectral region from 2 to 40 /im and are avail- 
able for Kl, K2 and K5 giants. Thus, after the spectral 
type was determined based on the K-hand spectrum, the 
Decin & Eriksson model closest to that spectral type was 
chosen. The model was reddened using the Indebetouw 
et al. (2005) extinction law as described above. After 
convolving the reddened model flux with the appropri- 
ate photometry point, the model was normalized to the 
data. In most cases, the reddened model was normalized 
to the Kg band, except for Q22-3 and Q21-4, where the 
IRAG4 and IRAG3 photometry, respectively, was used 
to normalize the reddened model. The reddened mod- 
els provide reasonable continua over the entire spectral 
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Figure 2. [top] H — Ks vs. J — K and [bottom] and J — K 
vs. Ks — [4.5] for the IC 5146 field stars (diamonds). In the top 
panel, smooth lines denote the dwarf (upper curve) and giant (lower 
curve) branches of normal unreddened stars. In the bottom panel, 
the curve for normal giant stars is shown. The reddening vectors 
are equivalent to Ay ~ 5 mag. 

energy distribution (Fig. [1|). The optical depth spectra 
were computed using r = — In F^/Fcontinuum, and are 
presented in Fig. |4l 

3. ICE AND DUST ABSORPTION FEATURES 

Several dust and ice absorption features are apparent 
in the optical depth spectra (Figure SJ: the 3.0 /xm II2O- 
ice and associated wing, the 6.0 /im and 6.8 /im ice/dust 
features, and the 9.7 /im silicate feature. Figure [7] shows 
the CO-ice absorption profllc for Q21-1. 

3.1. H20-ice and Long-wavelength Wing 

Water-ice is the primary ice-mantle constituent in 
dense clouds and its presence is evident through detec- 
tion of its fundamental 0-H stretching-mode absorption 
feature at 3.0 /xm and bending and libration modes at 
6.0 and 13 /xm, respectively. In laboratory ice analogs, 
the peak depth of the 6.0 /xm ice absorption feature is 
7 to 9 times weaker than the 3.0 /im absorption, and in 
the astronomical spectra absorption at 6.0 /im is a blend 
of H20-ice with other ice components (see below). The 
libration band is also weak and blended on the shoul- 
der of the broad 9.7 /im silicate absorption feature. The 
presence of alcohols and other compounds containing OH 
groups in the ice can contribute to the 3.0 /im absorption 
band and should be taken into account for those lines of 
sight that have GII3OII abundances greater than 20% 
relative to H2O. Since the methanol abundance for the 
objects studied here is on the order of a few percent, the 
contribution of alcohols to the 3.0 /im absorption can be 
neglected. Since the 3.0 /im band doesn't usually suf- 
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Table 2 

Peak Optical Depths of Various Ice and Dust Components 



Source ID 


T3.05 


''"3.47 


f"6.0 


comp. 1 


comp. 2 


comp. 3 


comp. 4 


comp. 5 




Q21-1 


1.45 


0.065 


0.19 


0.04 


0.04 


0.07 


0.06 


0.00 


0.602 ± 0.002 


Q21-2 


0.57 


0.030 


0.09: 


0.10: 


0.04: 


0.07: 


0.06 


0.00 


0.534 ±0.110 


Q21-3 






0.04: 












0.444 ± 0.026 


Q21-4 


0.31 


< 0.02 


0.01: 












0.281 ±0.082 


Q21-5 


0.30 


< 0.01 


0.04 


0.02 


0.02 


0.02 


0.02 


0.00 


0.259 ± 0.006 


Q21-6 


1.17 


0.055 


0.17 


0.05 


0.04 


0.08 


0.07 


0.00 


0.689 ± 0.008 


Q22-1 


0.82 


0.048 


0.11 


0.03 


0.04 


0.04 


0.04 


0.00 


0.509 ± 0.006 


Q22-3 






< 0.01 












< 0.15 


Q23-1 


0.24 


< 0.04 


0.06 


0.03 


0.03 


0.03 


0.03 


0.00 


0.297 ±0.017 


Q23-2 






0.05 


0.03 


0.03 


0.05 


0.05 


0.00 


0.379 ± 0.026 
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Figure 3a. Spectral energy distributions for the IC 5146 field 
stars with available IRTF-SpeX data. 2MASS photometry (di- 
amonds), SpeX spectra of IC 5146 field star (black line), SpeX 
standard star (green line) 

fer from severe issues of blending and is the strongest of 
the H20-ice absorption features, we use it to deduce the 
H20-ice column density as well as the ice "threshold" 
extinction. 

Fig. H] shows the laboratory spectrum of 10 K H2O- 
ice (Hudgins et al. 1993) overlaid on the optical depth 
spectra. Compared to the laboratory ice analog, the ob- 
served 3 /im ice profile has "excess" absorption longward 
of 3.05 /im and, in some cases, weak excess absorption 
at about 2.95 /im. Smith et al. (1993) noticed the latter 
feature in Taurus field stars on later spectral types, and 
therefore attributed it to a photospheric or circumstellar 
feature. However, for the IC 5146 stars, the spectra with 
the most apparent excess at 2.96 //m are late G and early 
K giants, so a stellar origin is less likely. Also present 
in these spectra is "excess" absorption in the 3.2 to 3.6 
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Figure 3a. continued 
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Figure 3a. continued 

/im region, known as the (enigmatic) "long-wavelength 
wing" (Baratta & StrazzuUa 1990; Smith et al. 1993). 
Ammonia (NH3) ice related species such as ammonium- 
hydrates may account for some of this absorption as well 
as the excess at 2.96 /im (e.g. Chiar et al. 2000; Gibb et 
al. 2001). Some of the wing absorption may be accounted 
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Figure 3b. IC 5146 field stars witli no available SpeX data. 

for by scattering in a model where spherical grains are 
coated with thick mantles or only the largest grains have 
mantles (Smith et al. 1993). Discrete features in the 3.2 
to 3.6 /im region contribute to the absorption: a fea- 
ture at 3.47 /im is attributed to hydrocarbon-containing 
ices (Chiar et al. 1996; Brooke et al. 1996, 1999) and 
CHaOH-ice exhibits a narrow absorption feature at 3.54 
fiui (AUamandola et al. 1992). However, in quiescent 
dense clouds the abundance of CH3OH relative to H2O- 
ice is less than 5%, so it does not contribute significantly 
(Chiar et al. 1996). Boogert et al. (2011) find the CH3OH 
abundance to be ~ 7% for the isolated dense cores probed 
by background stars. Thus, much of the "wing" absorp- 
tion remains a mystery. 

Four of our sources have a detectable 3.47 feature. 
The optical depth, T3.47 was calculated using a local third 
degree polynomial to determine the continuum across the 
3.2 to 3.8 /im region, avoiding the 3.30 to 3.33 /xm region 
(where the telluric cancellation is poor) and the 3.4 to 
3.6 /im region where the ice absorption is expected to be 
present. We determined T3.47 using the average feature 
width (AA — 0.105 ± 0.004/im) and central wavelength 
(Aq — 3.469 ± 0.002 /im) determined by Brooke et al. 
(1999). The results are given in Table 2. In addition, for 
Q21-1 and Q21-6, absorption in excess of the typical 3.47 
/im feature is present. Fig. [5] shows that this absorption 
is consistent with the CH3OH amounts allowed by the 
observed 9.75 /im feature in the spectra of these same 
sources. 

In spite of the excess wing absorption, the 3.0 /im ab- 
sorption feature itself is well-suited for calculating the 
Il20-ice column density. We estimate the H2 0-ice col- 
umn by scaling a Gaussian function with Aq = 3.04 /im 
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Figure 4. Optical depth profiles for tfie ice and dust features 
based on the reddened photosphere model matches as shown in 
the previous figure. Laboratory data (red lines) for H2O ice at 
lOK are overlaid on the observational data (black lines). Sources 
are (from top to bottom) Q21-1, Q21-6, Q22-1. 

with FWHM = 0.22 /im. On average, the column den- 
sities determined in this manner are about 10% larger 
than those determined using the 10 K H20-ice labora- 
tory spectrum. The resulting column densities are listed 
in Table [31 The laboratory ice analog spectrum is also 
useful for assessing the contribution of pure H20-ice to 
the 6.0 and 13.0 /im spectral regions. In order to limit the 
contribution of ice to those spectral regions, we fitted the 
laboratory spectrum to the trough of the 3.0 absorp- 
tion and calculated a residual across the entire Spitzer 
IRS spectrum. This calculated residual was the n use d to 
analyze the 9.7 /im silicate absorption profiles ( ii3.3p and 
the absorption features in the 5 to 8 /im region ( ij3.4p . 

3.1.1. H20-ice Threshold 

The ice threshold is the Ay value below which ice ab- 
sorption is not detected and therefore icy mantles are 
not present (Whittet 2003). The threshold Ay is gener- 
ally higher for dense clouds that have active star forma- 
tion and are forming intermediate- to high-mass stars, 
compared to quiescent clouds that form low-mass stars 
at a slow rate. Volatile ice species, like CO-ice, whose 
presence is more sensitive to local temperature and den- 
sity conditions tend to have higher threshold extinction 
(Whittet et al. 1989; Chiar et al. 1995). The Taurus 
cloud, generally used as a model for the measurement of 
ice thresholds in pristine environments since it has the 
greatest number of ice absorption measurements toward 
background stars, exhibits a threshold of Ay = 3.2 ±0.1 
mag (Whittet et al. 2001 and references therein). For 
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Figure 4. Profiles for (from top to bottom) Q21-2, Q21-3, Q23-2. 
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Figure 4. Profiles for (from top to bottom) Q21-4, Q21-5, Q23-1. 
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Figure 5. Absorption in tfie long-wavelength H20-ice wing [left 
panel] and 8.5 to 10 /im region after fitt ing local con tinua to remove 
H2 0-ice and silicate absorption (see ^3.11 and 13.51 for details) for 
Q21-1 [top] and Q21-6 [bottom]. The left panel shows a Gaussian 
representing the average 3.47 /im feature profile as determined by 
Brooke et al. 1999 (dashed line), the weak CH3OH ice feature 
(solid line) from laboratory data (Gerakines et al. 1995, 1996), 
and coaddition of those two features (deepest smooth line). The 
contribution of CH3OH at 3.54 fim is constrained by the strength 
of the maximum allowed CH3OH feature detected at 9.75 fim (left 
panel; laboratory data shown as smooth curve). 
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Figure 6. Optical depth of the 3.0 fim H20-ice feature versus vi- 
sual extinction for field stars behind the IC 5146 (this work, black 
symbols) and Taurus (Whittet et al. 2001, purple symbols) dark 
clouds. Ay values for the IC 5146 background stars are based on 
E{J — K) as determined from the redde ned p hot osph eric model fit 
and an appropriate extinction law (see i|2.1l and l2.2|l . Determina- 
tion of Ay for Taurus sources is discussed in detail in Whittet et 
al. (2001). A weighted least-squares fit to the IC 5146 data (black 
line) and Taurus (purple line) taking uncertainties in ra.o and Ay 
into account results in threshold extinctions of Ay = 4.02 mag 
and Ay = 3.19 mag, respectively. Accounting for the foreground 
extinction toward IC 5146 (~ 0.8 mag) brings the threshold Ay 
for IC 5146 and Taurus into agreement. 

the more active Serpens cloud, the H20-ice threshold is 
higher at Ay ~ 6 (Eiroa & Hodapp 1989), and for the p 
Ophiuchus cloud that is forming intermediate- to high- 
mass stars, the threshold is even higher at Ay ~ 13 mag 
(Tanaka et al. 1990). 

For the IC 5146 cloud, we calculate the H20-ice thresh- 
old extinction by means of a linear least-squares fit in- 
cluding the uncertainties (error bars) in both Ay and T3.0 



8 



Chiar et al. 




4.60 4.62 4.64 4.66 4.68 4.70 4.72 4.74 8 9 10 11 12 13 14 

Wavelength (micron) Wavelength (/^m) 



Figure 7. Optical depth spectrum in the region of the CO ice 
absorption feature for Q 21-1 (solid line). The (scaled) profile for 
the Taurus field star, Elias 16 is also shown for comparison (dashed 
line). See text for details. 

, resulting in a threshold of = 4.03±0.05 mag (Fig.lHl). 
For consistency, we repeat the fit to the Whittet et al. 
(2001) data while taking the uncertainties in both Ay 
and Ta.o into account, resulting in threshold extinction, 
Ay = 3.19 ± 0.07. As demonstrated by the correlation 
lines in Fig. [51 the slopes of the IC 5146 and Taurus cor- 
relation lines are similar: 0.072 ±0.002 and 0.068 ±0.003, 
respectively, suggesting that H2 0-ice mantle growth with 
increasing extinction is similar in both clouds. This is in 
line with calculations that show that once the critical 
threshold Ay is reached - i.e., the first mono-layer of ice 
is formed on the silicate or carbonaceous substrate - the 
local infrared radiation field has little effect on subse- 
quent ice layer growth due to strong H-bonding between 
neighboring H2 0-ice molecules (Williams et al. 1992). 
The apparent increase of the threshold Ay in IC 5146 is 
likely to be an effect of foreground extinction given the 
distance of 950 pc to IC 5146 (§1). In fact, Neckel & 
Klare (1980) estimate the foreground extinction in the 
direction of IC 5146 to be ^ 0.8 mag, comparable to the 
difference in the threshold extinction between IC 5146 
and Taurus. 

3.2. CO and CO2 Ice 

Toward field stars behind dense clouds, the more 
volatile CO-ice ranges in abundance between 10 and 40% 
relative to H20-ice (Chiar et al. 1995; Whittet et al. 
2009). We were able to acquire an M-band spectrum 
with IRTF-SpeX for one of the IC 5146 sources, Q21- 
1. Figure [7] shows the 4.67 /im CO ice absorption profile 
with the spectrum of the Taurus field star Elias 16 (Chiar 
et al. 1995) overlaid for comparison. Many of the sharp 
absorption peaks in the Q21-1 spectrum that dip below 
the Elias 16 spectrum are due to (unresolved) gas phase 
lines. Since the resolution is not high enough to model 
these lines, we use the relatively smooth Elias 16 spec- 
trum to estimate the CO ice column density. The Elias 
16 CO profile, and therefore also its column density, is 
scaled by a factor of 1.05 to give a CO ice column den- 
sity for Q21-1 equal to 6.83 x 10^^ cm'^ (Table©. This 
results in CO/H2O column density ratio equal to 0.27, in 
the same range found for other dense cloud lines of sight. 

The bending mode of CO2 occurs at 15.3 /zm and is 



Figure 8. Comparison of the silicate profiles for highest extinction 
sources Q21-1, Q-21-6, Q22-1 (dashed lines), the low-extinction 
source Q21-5 (dotted line) and the diffuse ISM as probed by WR 
98a (solid line; Chiar et al. 2006). The spectra are normalized at 
10.1 iim. Small differences are noticeable on the blue wing of the 
feature; these may be due to photospheric absorption in the field 
stars. For the high extinction sources, moderate excess absorption 
at the bottom of t he f eature (~ 9.8 fim) is due to CHaOH-ice 
absorption (see Fig. I12II . 

contained in the IRS short-high spectrum for Q21-1 pre- 
viously presented by Whittet et al. (2009). They found 
that the profile shape of the 15.3 /xm CO2 feature is un- 
varying in the dark clouds sampled (Taurus, IC 5146, 
Serpens) despite some scatter in the abundance rela- 
tive to Il20-ice which may be due differences in how 
the elemental oxygen is distributed in the main oxygen- 
containing ice species (H2O, CO, CO2). They find that 
the CO2 abundance relative to II2O is ^ 0.2 for Taurus 
and - 0.3 for IC 5146 (Q21-1). 

3.3. Silicates 

Differences in the profile shape of the interstellar sili- 
cate feature observed in dense clouds and the diffuse ISM 
have been noted previously. Generally, emissivity curves 
based on observations of the M supergiant fi Cephei 
(Russell et al. 1975) and Orion's Trapezium region (For- 
rest et al. 1975) are used to fit the observed absorption 
features in the diffuse ISM and dense clouds, respectively 
(e.g. Roche & Aitken 1984; Whittet et al. 1988; Bowey 
et al. 1998). The Trapezium silicate profile peaks at 9.6 
/zm compared to 9.9 ^m for the fi Cephei profile. Fur- 
thermore, van Breemen et al. (2010) find that the "dense 
cloud" silicate profiles differ substantially from the dif- 
fuse ISM profiles peaking at shorter wavelengths relative 
to the diffuse ISM profiles. To illustrate the robustness 
of these differences, we compare the IC 5146 field star 
profiles with the highest S/N to that of the diffuse ISM 
toward the Wolf-Rayet star, WR 98a (from Chiar & Tie- 
lens 2006) in Fig. [51 All profiles have been normalized 
to unity at 10.1 /im. The normalization wavelength was 
chosen to avoid regions that could be contaminated by 
photospheric or ice absorption. This comparison shows 
that there are small differences in the dense cloud spectra 
on the blue side; these differences may be due to residual 
photospheric absorption (like SiO) that is not fully ac- 
counted for in the models. In addition, there is "excess" 
absorption in the trough of the feature in the highly red- 
dened sources due to CII3OII- and NHa-ice absorption 
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Figure 9. The relationship between the optical depths of the 3.05 
and 6.0 fim features. Data points are shown for the observed fea- 
tures in IC 5146 (black points) and low-mass YSOs (green points, 
Boogert et al. 2008). The IC 5146 points are clustered in the lower 
left corner of the plot; the inset shows this region in more detail. 
The purple line is the ratio, rg.o/ra.os, for laboratory H20-ice at 
lOK (Hudgins et al. 1993). The black line is the least-squares fit 
to just the IC 5146 sources. 



(see tj3.5p . Otherwise, the silicate profiles show little dif- 
ference in the less- and more-reddened hnes of sight in 
this small sample. All the IC 5146 field star profiles are 
shifted blueward compared to the diffuse ISM profile. 

3.4. Absorption Features in the 5 to 8 fj,m Region 

Absorption features centered around 6.0 and 6.8 fim 
are detected in six out of ten lines of sight in IC 5146 (see 
Table 2). These features have been previously detected 
in the spectra of high-mass (Schutte et al. 1996; Keane 
et al. 2001) and low-mass young stellar objects (YSOs; 
Boogert et al. 2008), dense clouds as probed by a few 
field stars (Knez et al. 2005) and dense cores as probed 
by field stars (Boogert et al. 2011). 

It has been shown previously that the absorption fea- 
ture at 6.0 fj,m. cannot be fully accounted for by H2 0-ice 
(e.g., Tielens et al. 1984; Keane et al. 2001; d'Hendecourt 
et al. 1996). In fact, if the 6.0 /im absorption were fitted 
with Il20-ice, the O-H stretching mode at 3.0 is over- 
estimated by up to a factor of 3 for high-mass young stel- 
lar objects (Schutte et al. 1996; Gibb et al. 2000; Keane 
et al. 2001). While CH3OH also has a strong band at 
3.07 /im (e.g., Hudgins et al. 1993) which can skew the 
apparent 3.0/6.0 ratio, this is unlikely to be a problem 
in quiescent lines of sight since the CH3OH abundance, 
determined by analysis of absorption at 3 .54 and 9.7 /im, 
is less than 5% of the H20-ice (see §3.51 and also Chiar 
et al. 1996). Other proposed contributors to the absorp- 
tion at 6.0 /im are formic acid (HCOOH), formaldehyde 
(II2CO), and the C-C stretch in aromatics. Combining 
our small field star sample with the larger sample by 
Boogert et al. (2008; 41 low-mass YSOs, 8 high-mass 
YSOs and 2 field starsQ 

we examine the relationship 
between the 3.0 /im H20-ice and the 6.0 /xm absorption 

^'^ For sources with extremely deep 3.0 ^m bands (rs o > 5), 
where the depths cannot be measured due to lack of detectable 
flux, Boogert et al. (2008) estimate T3.0 using the H2O libration 
mode at 13 ^m. This is the cause of the larger error bars on those 
data points. 
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Figure 10. Residual absorption in the 5.4-7.6 /im region after 
subtraction of pure H20-icc. [left] A local straight-line continuum 
(dashed line) is fitted across this spectral region (see text for de- 
tails), [right] Continuum-subtracted spectra. Thin lines represent 
components CI, C2, C3, C4 as described in the text. The C5 com- 
ponent is not present in these sources. Thick lines represent the 
sum of the four components. 

bands. Figure [HI shows graphically, that the observed 
6 /xm absorption cannot be fully accounted for by 10 K 
Il20-ice. The shallower line shows the Tg. 0/73. 05 ratio for 
the 10 K laboratory H20-ice analog from Hudgins et al. 
(1993). The steeper line is the correlation line for the 
IC 5146 field stars. Significantly, nearly all of the data 
points, including those of the field stars where little to no 
(energetic or thermal) processing of the ices is expected, 
lie above the laboratory ice line. At larger values of T3.05 
(which roughly scales with Ay), the scatter between the 
data points increases, indicative of increased alteration 
of the ices in the most heavily obscured YSOs. 

Boogert et al. (2008) recently decomposed the 6.0 and 
6.85 /im features into 5 components in addition to a con- 
tribution from H20-ice using spectra from a large sam- 
ple of low-mass and high-mass YSOs. They decompose 
the 6.0 iiui band into two components referred to as CI 
(A = 5.84 /im) and C2 (A = 6.18 /im). The CI com- 
ponent is explained mostly by solid HCOOH and H2CO 
with abundances up to 6% relative to H2O (Boogert et 
al. 2008). The C2 component is also due to a blend 
of several species, with NH3 being the most dominant. 
JVIonomers, dimers and small multimers of H2O mixed 
with CO2 as well as contributions from salts are also pos- 
sible (Boogert et al. 2008). The 6.85 /im absorption band 
shows great variation from source to source for the high- 
mass and low-mass YSOs that have been studied (Keane 
et al. 2001; Boogert et al. 2008). This band is also de- 
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Figure 11. The correlation of the CI, C2, C3 and C4 absorption 
components with H2 0-ice column density. The optical depths of 
the CI through C4 components are normalized to N(H20). Purple 
square are the YSO data from Boogert et al. (2008) and solid circles 
are date for field stars from IC 5146 (this work), Ellas 16 (Taurus, 
Boogert et al. 2008), EC118 (Serpens, Boogert et al. 2008). 

composed into a short wavelength (C3, A — 6.755 /.tm) 
and long wavelength (C4, A — 6.943 /im) component by 
Boogert et al. (2008). Both components are attributed 
to the ammonium ion, which would have to be produced 
by low-temperature reactions since both components are 
also detected in quiescent cloud spectra. Underlying the 
CI through C4 components is a fifth broad component 
(C5) whose origin may be related to processing of ices in 
the vicinity of the embedded YSOs (Boogert et al. 2008). 

We use the component definitions as determined by 
Boogert et al. (2008) to investigate their presence in 
the field star spectra. We show the component decom- 
position for the four sources with the highest S/N in 
Fig. [TUl To be consistent with the determination of the 
CI through C5 components in the studies by Boogert et 
al. (2008, 2011), we fit a local continuum across the 5 to 
8 /zm region. The shape of the 5 to 8 /im region (where 
the absorption bands of photospheric H2O and SiO ap- 
pear) of the observed spectra is not always well-fitted by 
the computational models, so a straight-line continuum 
is fitted across this region. This local continuum, fitted 
using points around 5.5, 6.5 and/or 7.5 /im, removes the 
apparent offset of the optical depth spectra (see the left 
panels in Fig. ITOl) . Following the subtraction of the local 
continuum, we decompose the 6.0 and 6.8 /zm features in 
their components (see the right panels in Fig. [TU)) . The 
C5 component is not apparent in any of our spectra, in 
line with the suggestion that this component is due to 
the presence of processed ices in the vicinity of embed- 
ded objects (Boogert et al. 2008). We investigate the 
correlation of the CI through C4 components with the 
H20-ice column density. Along with the CI through C4 
measurements for the IC 5146 field stars in this paper, 
we also include the two field stars (Elias 16, EC118) and 
the YSOs from Boogert et al. (2008). For consistency 
with the Boogert et al. studies (2008, 2011), we normal- 
ize the component depth (e.g., rci) to N(H2 0) and plot 
that against N(H2 0) to trace the component growth rel- 
ative to the ice column (Fig. [TT|). The field stars (solid 
circles) seem to follow the trends previously noted for 
the YSOs. We find that C2 shows a tight relationship of 



decreasing strength with increasing N(H2 0) suggesting 
that this component is less volatile than H20-ice itself. 
CI and C4 show more scatter,. It is interesting to note 
that the CI and C2 components seem to be of equal 
depth in each field star spectrum; this does not appear 
to be the case for the YSO spectra studied by Boogert 
et al. (2008) and is perhaps indicative of the sensitivity 
of one or both of these components to the UV or tem- 
perature processing that takes place in the environment 
around YSOs. The (normalized) C3 component shows a 
flat relationship with N(H20), linking it as a non-volatile 
ice component that it not sensitive to varying tempera- 
ture or energetic conditions of YSO environments. 

3.5. Weak Ice Features in the 8.5 to 10 fim region 

The 8.5 to 10 /im region contains the N-H vibration 
inversion mode of NH3 (d'Hendecourt & Allamandola 
1986), the C-0 stretching vibration of CH3OH (Sandford 
& Allamandola 1993) as well as the 0-0 symmetric and 
asymmetric stretching mode of O3 (Bennett & Kaiser 
2005). Substructure at these wavelengths appears at the 
bottom of the silicate feature for the most highly ex- 
tincted sources, Q21-1, Q21-6 and possibly Q22-1. We 
also examine the high S/N spectrum of Q21-5. For this 
source, which has relatively lower extinction (see Ta- 
ble [T|), we do not expect to see the weak ice absorptions 
in the 8.5 to 10 /im region. So, we use this source as a 
check that the substructure detected in the high extinc- 
tion sources is not due to an artifact in the data reduction 
process. For all four sources, we fitted and subtracted a 
local polynomial continuum in the 8.4-10.5 /xm region of 
the optical depth spectra in order to subtract the contri- 
bution of the broad silicate absorption. The continuum 
fits and resulting residuals are compared in Fig. [12] in 
the left and right panels, respectively. In the spectra of 
Q21-1, Q21-6, Q22-1, an absorption feature appears at 
9.0 /im with FWHM of about 24 cm~^. A second ab- 
sorption feature at 9.7 /im, with FWHM of ^ 38 cm~^ 
is also present (tentatively in Q 22-1). We discuss the 
likely identifications of there features below. 

We attribute the weak absorption at ~ 9.7 /im, present 
in the spectra of Q 21-1 and Q 21-6 (and tentatively 
Q22-1), to the C-0 stretching vibration in CHaOH-ice, 
that occurs at about 9.75 /im (Sandford & Allaman- 
dola 1993). As shown in Fig. [3 the CH3OH identifi- 
cation is supported by residual absorption in the long- 
wavelength wing of the 3 /im H2 0-ice feature that is 
not accounted for by the typical 3.47 /im "hydrocarbon 
ice" feature. Assuming a FWHM equal to 29 cm~^ and 
intrinsic strength of A = 1.8 x 10"^^ cm molecule"^ 
(d'Hendecourt & Allamandola 1986), we calculate the 
column density of CH3OH for the three lines of sight 
(Table El). The CH3OH/H2O ratio is in the range of 1-2 
% for these sources. This is well within the previous lim- 
iting value (< 5%) determined for the Taurus dark cloud 
(Chiar et al. 1996). 

The N-H inversion mode occurs near 9.1 /im, and its 
central wavelength is highly sensitive to matrix neigh- 
bors. Embedded in a polar matrix, the feature peaks 
at 8.97 /im, and in an apolar matrix, the feature peaks 
at 9.57 /im (Lacy et al. 1998). The observed 9.0 /im 
feature falls in the range of the N-H umbrella mode, al- 
though it is narrower than the available NH3 ice analogs 
which have widths equal to ~ 50 cm~^ (d'Hendecourt 
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Figure 12. [left] Silicate profiles for the sources (from top to 
bottom) Q21-1, Q21-6, Q22-1, and Q21-5. Light line is the local 
polynomial continuum used to subtract off the silicate absorption 
to extract the weak ice features, [right] Continuum-subtracted op- 
tical depth spectra in the region of the silicate absorption. Sources 
are the same as in the left panel. The approximate wavelength po- 
sitions of NH3 and CH3OH are shown. Note that the peak wave- 
length of NH3 is extremely sensitive to temperature and neighbor- 
ing molecules in the ice. 



& AUamandola 1986; Lacy et al. 1998). Given the dif- 
ficulty in extracting the ice feature from the bottom of 
the broad silicate absorption, it is possible that we have 
underestimated the width by a factor of 2. Taking this 
uncertainty into account and using an intrinsic strength 
of A = 1.7 X 10^^^ cm molecule^^ for the 1^2 umbrella 
mode of NH3 (d'Hcndccourt & AUamandola 1986), we 
compute column densities for Q 21-1, Q 21-6 and Q 22- 
1. These are given in Table [3l For all three sources, the 
column of NH3 relative to H20-ice ranges from 2 to 5%. 

Ozone (O3) ice has absorption features at 9.06 fim 
and 9.64 fun due to the 0-0 symmetric and asymmetric 
stretch (Bennett & Kaiser 2005). The 0-0 asymmetric 
stretch is ^ 80 times stronger than that of the corre- 
sponding symmetric stretch (Bennett & Kaiser 2005). 
Given the depth of the observed absorption centered at 
9.7 /im, 03-ice would contribute a negligible amount to 
the optical depth at 9.0 /ini. Thus, it is unlikely that O3 
contributes significantly to the ices in these lines of sight. 

4. EXTINCTION 

In the diffuse ISM, the extinction due to silicates (i.e., 
Tg,r) shows a tight linear correlation with Ay (Whittet 
2003 and references therein). Whittet et al. (1988) and 
Chiar et al. (2008) showed that, in dense clouds, the opti- 
cal depth of the 9.7 ^m silicate feature (rgiiicato) is weaker 
per unit Ay compared to the diffuse ISM. The Chiar et 
al. (2007) study analyzed 30 spectra toward six differ- 
ent dense clouds, including IC 5146. For that sample, 
r9.7 was computed using a low-order polynomial fitted 
over the 5 to 15 /im region and an estimate oi E{J — K) 
based on a likely range of spectral types for the back- 
ground stars. The T9.7 and E{J — K) estimates thus 
vary slightly from the more thorough analysis presented 
here; most values are within 20% of those determined 
previously. A recent reanalysis of the Chiar et al. (2007) 
sample reports that the total optical depth at 9.7 /im 
(T9.7) correlates with the extinction in the K-band, Ak 
(McClure 2009), and the data points lie along the same 
correlation line as for the diffuse ISM. However, there are 
some caveats to this illusory agreement in correlations. 
First, the total r at 9.7 /im includes the underlying "con- 
tinuum" extinction, the extinction due to the silicates 
themselves, and a small contribution from ices (for the 
high extinction sources) so a correlation between T9.7 and 
Ak does not speak directly to the relation between the 
silicate and K-band extinction. Second, for the diffuse 
ISM studies, the optical depth due to the silicates alone 
(i.e., Tsiiicatc) is shown to correlate with the visual ex- 
tinction (Roche & Aitken 1984; Rieke & Lebofsky 1985; 
Bowey et al. 2004; Chiar & Tielens 2006). In the diffuse 
ISM studies, continuum extinction due to the ISM dust 
and the dust shells of the Wolf-Rayet stars is accounted 
for by a power-law extinction curve or 1/A emissivity 
curve, respectively. 

We illustrate the relationship between the extinction 
quantities in Fig. [T31 Optical depth is related to the ex- 
tinction by T\ = — ln(10~"^^/^'^). The upper frame in 
Fig. [131 plots T9.7, the peak sihcate optical depth exclud- 
ing the continuum extinction, versus the E{J — K) color 
excess for the IC 5146 sources and diffuse ISM sources 
(from Whittet 2003 and references therein). As discussed 
previously by Chiar et al. (2007), the diffuse ISM sources 
follow a tighter and steeper correlation compared to the 
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Figure 13. [top] The optical depth of the siUcate feature at 9.7 
/im versus the E{J — K) color excess, [bottom] The total optical 
depth at 9.7 /^m (including continuum extinction) vs E{J — K). 
Green points are for the diffuse ISM (Whittet 2003 and references 
therein). In the bottom frames, the diffuse ISM r values (green 
circles) do not include continuum extinction (see text for details). 

IC 5146 sources. In the lower panel of Fig. [T31 we fol- 
low McClure's method and plot the total optical depth 
at 9.7 /im (including continuum and silicate extinction) 
for the IC 5146 sources and the diffuse ISIVI sources. The 
continuum extinction at 9.7 /xm is calculated using the 
Indebetouw et al. (2005) extinction law [Ax/Ak — 0.43 
at 8 /im), assuming A^j = Ag, where Ak is determined 
from the continuum fits shown in Fig. [T] The correlation 
lines are least-squares fits to the data and show that in- 
cluding continuum extinction shifts both the diffuse ISIVI 
and dense cloud points up in the T9.7 vs. E{J — K) plot. 

5. SUMMARY AND DISCUSSION 

Most comparisons of the relative abundances, temper- 
ature, and evolution of ices in star-formation regions are 
made with respect to the Taurus dark cloud which is 
used throughout the literature as the prototypical pris- 
tine environment for the study of dust and ices. This 
is because the fortuitous placement of bright field giants 
behind the Taurus cloud have made spectroscopic studies 
practical. In this paper, we present a spectroscopic study 
of the quiescent IC 5146 dark cloud which was made pos- 
sible by the sensitive IRS instrument on Spitzer and the 
SpeX instrument on the IRTF. We have shown that the 
quiescent IC 5146 dark cloud shares many of the charac- 
teristics of the Taurus cloud; we summarize our results 
here. 

IRTF-SpeX 2 /im spectra were used to classify the IC 
5146 background stars with stunning accuracy, thereby 



Table 3 

Column Densities of Various Ice and Dust Components 



Source ID 


N(H20)=' 


N(CO) N(NH3) 
xlO^'' cm-2 


N(CH30H) 


Q21-1 


25.2 


6.83 0.6 - 1.1 


0.24 


Q21-2 


9.84 






Q21-3 








Q21-4 


5.18 






Q21-5 


5.18 






Q21-6 


20.0 


0.4-0.8 


0.32 


Q22-1 


14.2 


0.3-0.7 


0.16: 


Q22-3 








Q23-1 


4.32 






Q23-2 









No SpeX data exist for Q21-3, Q22-3, Q23-2. 
^ colon denotes uncertain value 

allowing continua to be fitted over the full 2 to 20 /xm 
range, resulting in reliable determinations of the extinc- 
tion, ice and silicate absorption profiles. We show that, 
after taking foreground extinction into account, the II2O- 
ice extinction threshold in the IC 5146 dark cloud is 
nearly equivalent to that in Taurus, showing the uni- 
versality of the onset of ice mantle growth in regions far 
from internal sources of radiation. The growth of the 
icy mantles as measured by the increase of ra.o with re- 
spect to Ay is similar in both clouds as shown by the 
comparable slopes of the regression lines on a t^xi~-^v 
plot. 

We have compared ice abundances in the IC 5146 cloud 
to Taurus and low- and high- mass YSOs (Table 4). Pre- 
vious estimates of the abundances of CII3OII and NII3 in 
dense clouds were based on limits only (Chiar et al. 1996; 
Gibb et al. 2001). However, recent laboratory (Bottinelli 
et al. 2010) and computational simulations (Cuppen et 
al. 2009) suggest that the hydrogenation process that 
forms CH3OII should be efficient in cold environments. 
The recent observational study by Boogert et al. (2011) 
of ices in isolated dense cores finds CII3OII abundances 
in the 5 to 12% range, in line with the idea of efficient 
production of CH3OH in the coldest environments. Here 
we have shown that CH3OH and NH3 are present in the 
quiescent cloud medium at the 2 to 5% level relative to 
H20-ice. In general, CH3OH ice is more abundant in 
YSO environments compared to quiescent dark clouds. 
The NH3 abundances for the dense cloud and YSO en- 
vironments shows wide variation, perhaps reflecting the 
difficulty of the spectroscopic measurement rather than 
the sensitivity of NH3-ice itself to cloud conditions. In 
either case, NH3-ice is an important mantle constituent 
in all dense cloud environments, and therefore an im- 
portant reservoir of elemental N in these regions. The 
CO-ice abundance appears to be similar in the IC 5146 
and Taurus clouds, while the CO2 abundance is elevated 
in IC 5146 compared to Taurus. Whittet et al. (2009) 
attribute this difference to variations in gas-phase chem- 
istry, but point out that the amount of oxygen depleted 
into the ices is consistent in the two clouds. As has been 
previously noted (Kerr et al. 1993; Thi et al. 2006), the 
YSOs show a large variation in CO abundance, reflect- 
ing the sensitivity of CO to the environmental conditions 
(temperature, UV), while the CO2 abundance is remark- 
ably consistent (with the exception of the IC 5146 mea- 
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Table 4 

Relative Ice Abundances in Different Environments 



Ice 



IC 5146 



Taurus 



low-mass YSOs 



high-mass YSOs 



H2O 

CO 

CO2 

CH3OH 

NH3 

OCN- 



1 

0.27 
0.35 
~ 0.02 
~ 0.05 
< 0.01 



1 

0.32 ±0.12 
0.19 

< 0.05 

< 0.10 

< 0.01 



1 

0.32 
0.19 ±0.13 
0.08 

< 0.27 
< 0.002 - 0.057 



1 

0.12 ±0.07 
0.18 ±0.07 

< 0.04 - 0.32 

< 0.02 - 0.17 
< 0.002 - 0.063 



surement!). 

Finally, our results firmly support the conclusions by 
Chiar et al. (2007) that the silicate optical depth, rg.r, 
grows more slowly per unit Ay in dense clouds relative 
to the diffuse ISM. The change in dust characteristics 
that cause the difference in the Tg,r-Av relation is still 
an open question. 
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